Ultrafast Quantum Gates in Circuit QED 
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We present a method to implement fast two-qubit gates valid for the ultrastrong coupling (USC) and deep 
strong coupling (DSC) regimes of light-matter interaction, considering state-of-the-art circuit quantum electro- 
dynamics (QED) technology. Our proposal includes a suitable qubit architecture and is based on a four-step 
sequential displacement of the intracavity field, operating at a time proportional to the inverse of the resonator 
frequency. Through ab initio calculations, we show that these quantum gates can be performed at subnanosecond 
time scales, while keeping a fidelity above 99%. 



Introduction. — With the advent of quantum information 
science [II, there have been enormous efforts in the design of 
devices with high level of quantum control and coherence 13. 
Circuit QED |I3}Q has become a leading technology for solid- 
state based quantum computation and its performance is ap- 
proaching that of trapped ions |6 | and all-optical implementa- 
tions 17 1 . Considerable progress has been made in recent cir- 
cuit QED experiments involving ultrastrong coupling 18] ID, 
two-qubit gate and algorithms |11-16|, three-qubit gate and 
entanglement llT7l - [T9l . Most of the proposed gates are based 
on slow dispersive interactions or faster resonant gates, and 
would require operation times of about tens of nanoseconds. 

To speed up gate operations, the latest circuit-QED tech- 
nology offers the USC regime of light-matter interactions |8 - 
UO 20, 21 1, where the coupling strength g is comparable to the 
resonator frequency ojr (0.1 < g/(^r ^ !)■ This should open 
the possibility to achieve fast gates operating at subnanosec- 
ond time scales |22 , 23| . In this sense, the design of these 
novel gates becomes a challenge as the rotating-wave approx- 
imation (RWA) breaks down and the complexity of the quan- 
tum Rabi Hamiltonian emerges [24, 25 [. Preliminary efforts 
have been done in this direction involving different configu- 
rations of superconducting circuits [26-28[. Likewise, in a 
recent contribution, it has been discussed the possibility of 
performing protected quantum computing l29l . 

In this Letter, we propose a realistic scheme to realize 
fast two-qubit controlled phase (CPHASE) gates between two 
newly designed flux qubits |30l , coupled galvanically to a 
single-mode transmission line resonator (Fig. [T]). Our pro- 
posal includes: (i) a CPHASE gate protocol operating at times 
proportional to the inverse of the resonator frequency; (ii) 
the design of the qubit-resonator system, allowing for high 
controllability on both the qubit transition frequency and the 
qubit-resonator coupling, in USC [8, 9[ and potentially the 
DSC regime [31 [ of light-matter interaction. Through ab ini- 
tio numerical analysis, we discuss the main features of this 
scheme in detail and show that the fidelity could reach 99.6%. 
This is an important step in the reduction of resources require- 
ment for fault-tolerant quantum computation [32[. 

Design of a versatile flux qubit. — The junction array is 
schematically depicted in Fig. [T] It consists of a six- 



Josephson-junctions configuration, each one denoted by a 
cross, coupled galvanically l8l|9l[33][ to a coplanar waveguide 
resonator. This is characterized by a superconducting phase 
difference Aip and the electromagnetic field supported. The 
upper loop stands for a three -junction flux qubit [30[, while 
the additional loops will allow a tunable quit-resonator cou- 
pling strength. We analyze this qubit design step-by-step. 
First, we describe the potential energy coming from the in- 
ductive terms, which is the dominant contribution. Second, 
we introduce the transmission line resonator inserted with a 
Josephson junction in the central line and introduce the super- 
conducting phase slip Atp across the junction shared with the 
resonator. Third, we add the capacitive terms which appear 
in the junctions and obtain the full Hamiltonian of the system 
using a standard procedure. Finally, we identify two levels 
in the degrees of freedom of the junction architecture, which 
will define our qubit, and obtain the effective Hamiltonian to 
describe the fast two-qubit gate. 

The potential energy due to the inductive terms is ob- 
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FIG. 1. (Color online) Circuit-QED configuration of a six- 
Josephson-junctions array coupled galvanically to a resonator (bot- 
tom line). The flux qubit is defined by three Josephson junctions in 
the upper loop threaded by an external flux <l>i . Two additional loops 
allow a tunable and switchable qubit-resonator coupling. 
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tained by adding up the corresponding Josephson potentials 
£{tfk) — —Ejk cos{(pk), where Ejk and (pk represent the 
Josephson energy and the superconducting phase across the 
k-th junction, respectively. We assume Eji = Ej2 = Ej, 
Ej3 = aEj and Ej4 = Ej^ = ctiEj. In addition, around 
each closed loop the total flux has to be a multiple of the flux 
quantum $o = h/2e, or expressed in terms of superconduct- 
ing phases, J^kfk = ^irfj + 27rn, where we defined the frus- 
tration parameter fj — $j7$o = 1, 2, 3). Using this quan- 
tization rule, the total potential energy U, containing both the 
qubit energy and the qubit-resonator interaction, reads 

= — [cos(y5i + cos(y32 + acos((^2 — + 2nfi) 

Ej 

+ 2a4(/3)cos(¥>2-<^i+2^/ + A^/.)], (1) 

where q;4(/3) = "4 cos(7r/3), / = /i + /2 + /3/2, and Aip 
stands for the phase slip shared by the resonator and /2-loop. 
Note that the junction at the central line introduces a boundary 
condition that modifies the mode structure of the resonator, 
but without altering the inductive potential The parame- 
ters a, Q!4 and fj can be optimized to find a suitable working 
point. 

We introduce now the Hamiltonian of the interrupted in- 
homogeneous transmission line resonator After mapping the 
nonuniform resonator into a sum of harmonic oscillators, the 
standard Hamiltonian described by the flux amplitude tpn and 
charge g„ as its conjugate momentum is ETI 

where C,. = Cr + Cjq is the modified resonator capacitance 
due to the presence of the sixth Josephson junction at the cen- 
tral line, see Fig.[T] The quantization procedure yields 

Hr = huJr, (a+a„ + 1/2) , (3) 

n 

with Ldn = kn/ ^/ LqCq (Cq, Lq the capacitance and induc- 
tance per unit length, respectively), and fc„ — {2Lo/Lj){l — 
ujf^/ujp) cot(fc„l), where ujp — 1/y/LjCj is the plasma fre- 
quency of the junction and 2 1 is the length of the central line. 

For simplicity in the presentation, assume that the junc- 
tions only interact with the first eigenmode of the resonator 
The modified superconducting phase slip Atjj is related to the 
single-mode resonator variables as 

A^:^ At/jiia + a^), (4) 

where Aipi = {Si/ipo){^/2cOrCrV^^ ■ Here, ipo = ^o/2tt 
is the reduced flux quantum, ujr is the frequency of the first 
eigenmode. Si = ui{x2) — ui{xi) corresponds to the differ- 
ence between the first-order spatial eigenmode, evaluated at 
points shared by the resonator and the /2-loop. 

The total Hamiltonian is obtained by including the kinetic 
terms of the qubit and the Hamiltonian of the resonator 

n = 4:AEc{nl + nl) + 8BE^nin2 + Hr 
C 

+ 2e—qn{ni- n2) + U{Lpi,ip2), (5) 



where Ec = e^/2Cj is the charging energy. A, B, C are func- 
tions of system parameters |34|, such as junctions size and 
phase slip magnitude. The degrees of freedom of the junction 
architecture are {tpi , (P2) and their conjugate momenta, which 
are the numbers {ni ,712) of Cooper pairs in the boxes. 

In order to define a qubit within the junction architecture, 
we diagonalize the term of the Hamiltonian containing only 
the junctions. The two lowest energy eigenstates are labeled 
as the eigenstates of cr^, and the two-dimensional subspace 
spanned by them describes the qubit. Furthermore, since 
AV'i ^ 1 in general, we expand the potential ([T]l up to the 
second order in Atjj. This gives rise to the first-order and 
second-order qubit-resonator inductively coupling. Finally, 
after projecting the interaction terms also into the qubit basis, 
the Hamiltonian reads 

v. = (Jz + hujra^a + Hint (6) 
with the effective interaction Hamiltonian 

Hint - 2i?ja4(/3) ^ (A^)" J2 c;r("'"4,/l,/2K, 
m— 1,2 fi—x,y,z 

(7) 

c™{a,a4, fi, f2) being the controllable magnitudes of the 
longitudinal and transverse coupling strengths for m-th order 
interaction. Here we have ignored the capacitive coupling, as 
it is orders of magnitude smaller than the inductively coupling. 

Numerical analysis. — We provide ab initio numerical ex- 
amples to show the functionality of our setup. Firstly, we 
study the properties of the inhomogeneous transmission line 
resonator obtained by inserting the sixth Josephson junction 
on the central line. The eigenfrequencies w„ and eigenmodes 
Un{x) can be found numerically and calculated with indepen- 
dence of the qubit Hamiltonian f2T|. In fact, considering a 
resonator of impedance Z ^ 50 ft, capacitance Cr ~ 850 fF, 
and a Josephson capacitance Cjg ~ 10 fF, we estimate the 
first mode of the resonator with frequency LOr /2tt ~ 7 GHz, 
which leads to a phase slip magnitude A-ipi — 0.1218. These 
values determine the qubit-resonator coupling strength. Sec- 
ondly, we perform the numerical study of the total Hamil- 
tonian (jSj), which shows that when external fluxes satisfy 
/2 + /3/2 = 0.5, both and the second-order coupling are 
negligible. This reduces the interaction Hamiltonian to 

"Hint = %(a + a^)(cz(Jz + Cx<Jx), (8) 

with the effective coupling strength g = 2Ej a^ifs) Aipi/h 
and Cz^x = cl z- 

In Fig. |2|a,b) we plot and as a function of the qubit 
junction size a, and the frustration parameter fi. They clearly 
show another characteristic of the setup, that is, the switching 
from transversal to longitudinal couplings depending on the 
external flux $1 iflOl . In particular, when selecting a qubit 
junction size a = 1.2 and /i = 0.505, we obtain a large 
(small) contribution of longitudinal (transversal) coupling — 
see Fig. ^c) where Cz and Cx are depicted for a parameter 
/a = and /g = 1. 
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FIG. 2. (Color online) (a,b) Couplings strength Cz and Cx as a func- 
tion of the qubit junction size a = Ejs/Eji, and the external frus- 
tration parameter /i — ^i/^o, where "l>o is the flux quantum, for a 
frustration /a — 0. (c) Coupling strengths for a qubit junction size 
a — 1.2, and for values fa — and fs — 1. (d) Qubit spectrum 
for values fg = {1,0,0.5}. In this simulation we have considered 
Ej/h = 221 GHz, a junction size 04 = 0.058. 

From the diagonalization of qubit Hamiltonian for different 
values of the frustration parameter /a = {0, 1, 0.5}, we esti- 
mate a qubit frequency ujg/2TT - {10.94, 11.25, 10.99} GHz 
[see Fig. |2|d)]. Furthermore, considering a junction size 
a4 = 0.058, the qubit-resonator coupling strength is g/uJr — 
{0.446, —0.446, 0}, reaching the USC regime. It is notewor- 
thy to mention that our setup allows us to turn on/off the cou- 
pling g as well as to change its sign, operations that may be 
carried out in times of the order of 0.1 ns, or even less |[35l[36l . 

The ratio g/ujr could change at will without affecting sig- 
nificantly the qubit properties required for the proposed pro- 
tocol. For instance, for a junction size a4 = 0.12, a phase 
slip magnitude Aijji = 0.0768 and a resonator frequency 
ujr/^TT ~ 8.01 GHz (obtained from a central Josephson ca- 
pacitance C,/6 ~ 17 fF), one obtains g/ujr = 0.509. Also, 
with other choices of parameters, we can move from the USC 
to DSC regime and even the ratio of transversal to longitudi- 
nal coupling could be tuned. These examples, and particularly 
the model introduced by Eq. ([8]l, will be the basis to develop 
protocols for fast two-qubit gates. 

Fast two-qubit gate. — The external tunability of the previ- 
ous circuit is now exploited to propose an fast protocol for 
a two-qubit gate ll22l . In Ref. |26| a scheme was studied, 
where a two-qubit CPHASE gate is produced by alternating 
between positive and negative values of the coupling strength 
g for each qubit. In our architecture, this can be done simply 
by changing the flux /a from to 1 ; however, this action will 
also increase the value of the undesired transversal coupling 
in Hamiltonian of Eq. (j6]l. For instance, in case of a junc- 
tion size = 0.12, we find g/ujr = {0.509,-0.509} and 



= {0.040, 0.135} for /3 = {0,1}. 

Hence, a more suitable protocol for the proposed architec- 
ture consists of: 

Step 1.- The coupling gi is maximized (f^^^ ~ 0), whereas 52 
made exactly zero by tuning /g ' — 0.5. The system 
evolves for a period ojr^i G (0, 7r/2]. 

Step 2.- The coupling 172 is maximized (/g = 0), whereas gi 
made exactly zero by tuning /g^^ = 0.5. The system 
evolves for a period ujrt2 = tt — uirti. 

Step 3.- Repeat Step 1. 

Step 4.- Repeat Step 2. 

We study first the ideal case, in which the transversal com- 
ponent of the coupling in Eq. (j6]l is negligible, which could be 
achieved by tuning the fluxes /i = 0.505 (or larger) of each 
qubit. Then, the two-qubit Hamiltonian is 

H^Y. ^^i*' + ^^ra'^a - ^ hg,{a + a^)cri'l (9) 

i i 

Under this Hamiltonian, the unitary evolution operator Ui cor- 
responding to each step is 

\UJr 

U2,4 = 2 t2^-rujrt2a'ajy y2 ^^^^^t^ _ ^( 

\U)r 

where V {ficJz) — exp ({/3a''' — fi*a}(7z^ is a controlled co- 
herent displacement of the field and ojg^^ stands for the 
value of the qubit frequency when gi is maximum, zero, re- 
spectively. Finally, using V{a)V{(i) = e'^"'^°'^">V{a + /?) 
and e-^^"^"I?(a)e^^"^'' = r»(ae-'*), the gate W = JJ^U^ is 

K — g-»("91t*l+"9U*2)o-i^' 

,1., t 4-,. t ^rrm 4isinuj,ti^cri^^cri^^ 
X e'"^ ''24.*l+"92t*2)'^z g - ^ . (11) 

This gate U is equivalent to a CPHASE quantum gate |[T|, 
up to local unitary operations, provided 4 smuJrtigig2/(^r — 
7r/4 (notice that this condition requires gi/uJr to be in the 
ultrastrong regime). In the case of having a junction size 
04 — 0.12, for which the coupling strength takes a value 
of gi/iOr — gil'-^r — 0.509 and uj^. ^ Itx y. 8.01 GHz, 
then ijjrt\ = 0.86 and the total gating time will be tgatc — 
27r/wr — 0.12 ns. This is much shorter than typical coher- 
ence times in these systems, which are around 1 ^s. 

Deviations from perfect fidelity are expected if one ac- 
counts for undesired transverse coupling in Eq. (|6]l. For an ini- 
tial state where both qubits are in state |+) = + \e))/\/2 
and the resonator in the vacuum, we can compute the fidelity 
of the state generated assuming that cx' = 0.040 for each 
qubit at /g — 0. The fidelity of this state, with reduced 
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density matrix p, as compared to the ideal \tlj) {tp\, for which 
transverse coupHng is neglected, amounts to = {ipl p\ip) > 
0.996. This result is unchanged even if we include up to the 
third cavity mode in our ab initio calculation. For the sake of 
simplicity, we have considered instantaneous changes in the 
value of the fluxes /a. In this sense, the scheme can be easily 
adapted to account for smooth time-dependent profiles in the 
value of the coupling strength of both qubits, provided the ad- 
equate interaction time and number of iterations, and that no 
overlap between the pulses occurs. Indeed, switching frequen- 
cies of about 10 — 80 GHz are already available t35i i36l. This 
should allow the experimental realization of a high-fidelity 
fast CPHASE gate with subnanosecond operation time. 

Discussion. — Given the proposed CPHASE gate, based on 
the tunable qubit-resonator coupling in USC, we may con- 
sider the following extensions, (i) Multi-qubit entanglement 
and gate operations, such as realization of three-qubit Toffoli 
gates in the USC regime, faster than other schemes working in 
the strong coupling regime lfT8l[T9l[37l . (ii) With the advan- 
tage of switchable coupling in both strength and orientation, 
we may think of generating Ising-type Hamiltonians for qubit 
arrays, (iii) By controlling the geometric-related flux values, 
we can increase the higher-order couplings and thus study the 
nonlinear dynamics in USC regime, (iv) The adjustable cou- 
pling also allows us to couple to slower measurement devices. 

Conclusion. — We have proposed a realistic scheme for im- 
plementing an fast two-qubit CPHASE gate in current circuit- 
QED technology. The gate may work at subnanosecond time 
scales with fidelity — 0.996. This proposal may lead to a 
significant improvement in the operating time with respect to 
standard circuit-QED scenarios, as well as microwave/optical 
cavity QED systems, together with the reduction of the large 
resource requirement for fault-tolerant quantum computing. 
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